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Abstract 

The degrees of freedom of any interacting quantum field theory are entangled 
in momentum space. Thus, in the vacuum state, the infrared degrees of freedom 
are described by a density matrix with an entanglement entropy. We derive a 
relation between this density matrix and the conventional Wilsonian effective 
action. We argue that the entanglement entropy of and mutual information be- 
tween subsets of field theoretic degrees of freedom at different momentum scales 
are natural observables in quantum field theory and demonstrate how to com- 
pute these in perturbation theory. The results may be understood heuristically 
based on the scale-dependence of the coupling strength and number of degrees 
of freedom. We measure the rate at which entanglement between degrees of free- 
dom declines as their scales separate and suggest that this decay is related to the 
property of decoupling in quantum field theory. 



1 Introduction 



A quintessential property that distinguishes quantum mechanics from classical me- 
chanics is the entanglement of otherwise distinct degrees of freedom. When certain 
degrees of freedom are entangled with the rest of a quantum system, it is not possible 
to describe them by a pure state. Rather, the most complete description of a sub- 
system A at a particular time is via the reduced density matrix obtained by tracing 
over the degrees of freedom in the complement, pa = tr^(|\l/)(\l/|), where is the 
state of the entire system. The entropy constructed from the reduced density matrix, 
S(pa) = — tr(p^ln(p^)), quantifies the amount of entanglement between A and its 
complement. The entanglement entropies corresponding to reduced density matrices 
for diverse subsets of degrees of freedom provide a rich characterization of the quantum 
state for systems with many degrees of freedom^] 

In physical systems, we typically only have access to a subset of the degrees of 
freedom, namely the low-energy or long-wavelength modes that are directly accessible 
to experiments. In an interacting theory, it will generally be true that these degrees 
of freedom are entangled with the inaccessible high-energy degrees of freedom. Thus, 
the long-wavelength modes will be described by a density matrix. A more familiar 
description of low-energy degrees of freedom is due to Wilson [2] - one carries out the 
complete path integral over the inaccessible degrees of freedom, arriving at an effective 
action capturing the dynamics of the remaining system. We show that the Wilsonian 
prescription is compatible with the description in terms of a density matrix: given 
a Wilsonian effective action we can canonically associate the corresponding density 
matrix via equation ([5]) below. 

For continuous physical systems described by interacting quantum field theories, 
understanding the variation with scale of the Wilsonian effective action Sw{p) pro- 
vides key insights into the nature of the quantum field theories, revealing a striking 
insensitivity of the low-energy physics to the details of the ultraviolet description. Cor- 
respondingly, it is natural to consider the variation with scale of the density matrix 
p(/i) for the degrees of freedom with momentum \p\ < p and the associated entangle- 
ment entropy S(p). To make our considerations concrete, we derive a formula for this 
low-energy entanglement entropy in perturbative quantum field theory, and apply it to 
scalar field theories with <p n potentials in various dimensions. The scale-dependence of 
the entropy S(p) in such theories can be understood in terms of the variation of the 
coupling and number of degrees of freedom with scale. 

To study entanglement between scales in greater detail, we can consider the entan- 
glement entropy associated with any subset of the allowed momenta, or the mutual 
information between any two subsets of the allowed momenta, for example between 
individual modes with momenta p and q. These measures characterize the extent of 
entanglement between specific scales in field theory, and we compute the rate at which 
this entanglement declines as the scales separate. This fall-off may give an alternative 
characterization of the property of decoupling in quantum field theories. In theories 

^or a basic review of density matrices, entanglement entropy, and related concepts, see for example 
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that do not enjoy the property of decoupling, e.g. noncommutative gauge theories [3] 
and theories of gravity, the entanglement between degrees of freedom at different scales 
may play an especially important role. 

Entanglement entropy in quantum field theory has been considered previously, but 
almost all previous work has focused on entanglement between degrees of freedom 
associated with spatial regions (e.g. [U [5] ). The notion of a density matrix for low- 
momentum modes or entanglement between different momentum modes has appeared 
earlier in the context of cosmology and condensed matter physics (e.g. [61 E]), but 
there is little overlap with the present work. 



2 The low-energy density matrix 

A quantum system with many degrees of freedom has a Hilbert space of the form 
T-L = "Hi <8 H-2 ® ■ ■ ■ ■ Given a subset of these degrees of freedom (A, with complement 
A), we can write Ti = Ha^Ha where %a is the tensor product of Hilbert spaces for the 
degrees of freedom in A. If p is the density matrix for the full system (which may be in 
a pure state), a reduced density matrix for A is defined by tracing over A: p A = tr^(p) 
(or, given components in a specific basis, p^ n = J2 N p m N,nN)- Expectation values of 
operators that act only on A can be computed as tr (p {Oa <8> 1)) = ^a(p A Oa)- If A is 
entangled with its complement, p A will have a finite entropy: S = — tr a(p A \ogp A ) > 0. 

In this construction, the Hilbert space can be be decomposed into a tensor product 
in any convenient way. For example, the Hilbert space for two identical oscillators could 
be decomposed either as a product of the Hilbert spaces for the individual oscillators, 
or as a product of the Hilbert spaces of even and odd normal modes. A reduced density 
matrix could be computed in either case - good choices of decomposition are dictated 
by the structure of the interactions and restrictions on which degrees of freedom are 
accessible to measurements. 

In quantum field theories, locality makes it natural to associate independent de- 
grees of freedom with disjoint spatial domains. Hence, given a spatial region A (and 
complement A), one can decompose the Hilbert space as H = Ha <8> Ha an d trace 
over A to derive the reduced density matrix of A. But since the Hamiltonians of free 
field theories are diagonalized by modes of fixed momentum, it is in many ways more 
natural to use the Fock space decomposition, % = where is the Hilbert space 

of modes of momentum While this decomposition is motivated by considering the 
case of free field theory, it applies equally well once we turn on interactions, and is 
indeed the standard setting for computations in perturbative quantum field theorylf] 

2 Here, it is clearest to define the field theory as a limit of a theory on finite volume so that the 
tensor product is over a discrete set of allowed momenta. For a general held theory, the factors would 
be labeled by held type and spin/polarization in addition to momentum. 

There is a formal sense in which turning on interactions takes one out of the original Hilbert 
space. However, by placing a cutoff at some energy scale much higher than any scale of interest, the 
Hilbert space structure of the interacting theory will be the same as that of the free theory, and a 
density matrix for low-momentum modes can be precisely defined. Furthermore, as we will see later, 
various observables related to the spectrum of the density matrix have a well-defined limit as we take 
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In free field theory, the vacuum state is a tensor product of the Fock space vacuum 
states for each independent field mode - there is no entanglement between the field 
modes at different momenta. But in an interacting theory, the full vacuum state will 
be a superposition of Fock basis states - hence the modes of different momenta will 
generally be entangled. In this case, the reduced density matrix corresponding to a 
subset of the degrees of freedom (A) will necessarily have a finite entropy, indicating 
that A is effectively in a mixed state if the rest of system is traced over. Now, one 
is most often interested in the physics of the "infrared" degrees of freedom that are 
accessible to experiment, i.e., the degrees of freedom with momenta below some scale \x. 
The present discussion shows that tracing over the ultraviolet, i.e. degrees of freedom 
with momenta above /i, should lead to an infrared effective description in terms of a 
low- energy density matrix corresponding to a mixed state with finite entropy. 



Relation between low-energy density matrix and low-energy effective action 

The standard way of studying the low-energy theory is through the Wilsonian effective 
action. How is this quantity related to the low-energy density matrix? To begin, 
consider a bare action Sa defined with a cutoff \p\ < A. Associated to this, we have a 
Hamiltonian H\, which will have some ground state and corresponding density 
matrix p A = |^ A ) (^°|. This density matrix can be written as a Euclidean path integral 
by taking the T — Y limit of the finite temperature density matrix 

1 1 r<i>(r=(S/2)=$ y 

(<Py\p T M = Uhle-^M = tf / e~ S ° , (1) 

Z Z J <t>(r=- 13 /2)=<P y 

where {<ft y } is a basis of field configurations indexed by y, (3 — 1/T is the inverse 
temperature, is the Euclidean action, and Z is the partition sum that normalizes 
the path integral. 

Given a subset of degrees of freedom A (with complement ^4) and the tensor product 
structure of the Hilbert space, we can split the parameter y which indexes the basis 
states as y = (a, a), and a pick a basis with <p y = (p a x (f> a . To define a reduced density 
matrix for A by tracing over A we write 



<M/3/2)=0 a 



In the last expression, the fields <j)^ are periodic in the range [— /3/2, /3/2], which is 
implied after substituting (JTJ) into the trace in the middle expression. 

Now define a conventional thermal effective action for the the subsystem A: 

e -S£^ A ) = / V(p A {T) e -SE{J>A^A) . (3) 

^-/3/2<r</3/2 



the cutoff to infinity. 
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In terms of this, the density matrix for A is 



-, r <t> A (0/2) = $a 1 r <P A (T = 0+)=4>a 

, T U \ — _ / T)rh p ~ s w^a) — _ / TtAA 



(Ml^a) = -= / VMr) e~^ ((pA) = - / VMr) e 

L J4>A{-P/2)=<t>a L J4>A(r=0-)=<p a 

(4) 

In the last expression we translated time to put the discontinuity in the integral at 
r = ± and the fields are taken to be periodic at r = ±/3/2. The reduced density 
matrix for A in the ground state |\E^) for the entire system is extracted by taking 
->• oo. 

We now specialize to the case where A is the subset of degrees of freedom with 
spatial momenta \p\ < p for any scale p which is lower than the ultraviolet cutoff A. 
The reduced density matrix p\ v \ <it obtained by tracing over the degrees of freedom with 
momenta in the range p < \p\ < A is thus given by: 



I /'</'|p|<m( t = 0+ )=^IpI<m 

(4>\ p \<Mp\<M\p\<») = y Vcf> ]p]<tl e- s ^<J . (5) 

/ ^b|<M^ = ~) = *bl<M 



where now (having taken (3 — > oo,) Syy is the standard Wilsonian effective action 
obtained by integrating out the degrees of freedom with large spatial momentaQ 

e -s w (4> M< ^ = / D0| p | >M (r) e- SE{(t> ^<^>^ . (6) 



Equation (jSJ) is our final result for the low-energy density matrix. In particular, if 
O is an observable built out of the low- momentum modes at r = 0, it follows from (jSJ) 
that 

tr(0p) = || [d^ w< JOe-^»w<M) , (7) 

so the standard calculation using the effective action is equivalent to a calculation 
using the density matrix. Of course, the full Wilsonian effective action contains more 
information than the density matrix associated with the vacuum state of the field 
theory. The former is a functional of time-dependent field configurations, while the 
latter depends only on a pair of time-independent field configurations. 

The description of low-energy degrees of freedom via a density matrix may seem 
unfamiliar and one may ask why we cannot simply associate a pure vacuum state to 
the low-energy degrees of freedom based on the effective action. The reason is that Sw 
will typically contain terms with higher time derivatives, there is no way to associate to 
Sw a Hamiltonian expressed exclusively in terms of the low-momentum variables 
and their conjugate momenta. Thus, there is no canonical way to associate a pure 
state of the low-momentum part of the Hilbert space to the full ground state of the 
theory. As we have seen, the object that can be canonically associated to a Wilsonian 
effective action for these low-momentum degrees of freedom is a density matrix. 



4 Note that while the path integral is Euclidean, we are integrating out all modes with spatial 
momenta \p\ > /i, regardless of frequency. 
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3 Measures of entanglement 



What observables quantify the amount of entanglement between the degrees of freedom 
in different ranges of momenta? In this section we begin by discussing such quantities in 
generality and conclude by constructing perturbative expressions for such observables 
in weakly coupled field theories. 

First, for any density matrix p, the von Neumann entropy 

S(p) = -tr(pln(p)) (8) 

measures the classical uncertainty associated with the mixed state described by p. 
When p represents a microcanonical or canonical ensemble, the von Neumann entropy 
gives the thermodynamic entropy. When p is the reduced density matrix describing a 
subsystem A of a quantum system that is in a pure state, S quantifies the entanglement 
between A and its complement (A). In this case the entanglement entropy of A is equal 
to that of A, a fact that follows from a stronger result that the spectrum of eigenvalues 
of pa matches the spectrum of eigenvalues of px- 

When the Hilbert space for the theory can be decomposed into a tensor product 
with three or more factors, the quantum entanglement and classical correlations be- 
tween pairs of these subsystems are jointly quantified by the mutual information. For 
instance, if the Hilbert space is of the form % = Ha ®T-Lb® He <8> • • ■ , the mutual 
information between A and B is defined as 

I (A, B) = S{A) + S(B) - S(A U B) . (9) 

where S(X) is the von Neumann entropy of the reduced density matrix of subsystem 
X. Mutual information is always greater than or equal to zero, with equality if and 
only if the density matrix for the AB subsystem is a tensor product of the reduced 
density matrices for subsystems A and B. In other words, mutual information is zero 
if and only if there is neither any entanglement nor any classical correlation between 
the two subsystems!! Mutual information provides an upper bound on all correlators 
between the two regions: for any bounded operators Oa and Ob, acting only on the 
subsystems A and B, we have [S] 

J(A m > ((OaOb}-(Oa}(Ob)) 2 (m) 
/(A B) ~ 2\0 A \*\0^ • (10) 
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If the Hilbert space consists of three factors % = %a ® Hb <E> He an d the complete 
system is in a pure state it follows from the definitions that 

I{A U B, C) = I (A, C) + I(B, C) . (11) 



5 The systems A and B are said to be entangled if and only if the density matrix for the AB 
subsystem cannot be written as PiP\ ® Pb ■ Separable density matrices of this form represent states 
which have no quantum entanglement, but may have classical correlations. The mutual information 
for such a state can be nonzero. 



5 



But if A, B, and C together comprise only a part of the system, another interesting 
observable is the tripartite information which quantifies the extent to which the mutual 
information between A U B and C is determined by the parirwise mutual informations 
I(A,B) and I(B,C): 

I (A, B, C) = I (A UB,C)- I(A, C) - I(B, C) . (12) 

In general, this quantity can be positive, negative, or zero. For a pure state of the full 
system, I(A, B, C) is symmetric between the four subsystems A,B,C, and AU B U C . 



3.1 Entanglement observables in perturbation theory 

For weakly coupled quantum field theories, we can use perturbation theory methods 
to calculate the entanglement observables described in the previous section. To begin, 
it is useful to derive a set of perturbative results that apply more broadly. 

Consider a general quantum system whose Hilbert space may be decomposed into 
a tensor product % = Ha ® Hb, and start with a Hamiltonian of the form 

H = H A (gl + l® H B . (13) 

Denote the energy eigenstates of Ha by \n) and the energy eigenstates of H B by \N), 
with energies E n and En respectively. Before adding interactions, the ground state is 

|0,0) = |0)<g>|0). (14) 

Now, perturb the Hamiltonian by an interaction XHab, where A is a small parameter. 
The perturbed ground state may be written (before normalization) as 

i^> = io, o) + A >> °) + Yl Biv i°' N )+ Yl N ^ ' ( 15 ) 

where A,B, and C are coefficients starting at order A that may be computed in pertur- 
bation theory. To normalize, we should multiply by I/A/" 5 , where M = 1 + \A n \ 2 + 
^2 \Bn\ 2 + ^2 \C n ,N\ 2 ■ Now, the density matrix corresponding to the subsystem A is 

1 ( l + \B\ 2 A^ + BC^ \ , . 

PA 1 + \A\ 2 + \B\ 2 + |C| 2 V A + CB* AA^ + CC^J' [ } 

where the elements of this matrix correspond to |0)(0|,|0)(n|,|m)(0|,|m)(n| terms re- 
spectively. By a symmetry transformation p — > MpM -1 , we can simplify the form 
to 

PA=( 1 ~W* c ° c ^+0(\*) (17) 

where we are using the fact that A, B, and C start at order A. (See below for why this 
is possible.) 
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Up to corrections of order A 3 , the eigenvalues of this matrix are A 2 aj and 1 — A 2 £ a ii 
where {a{\ (normalized to be of order A ) are the eigenvalues of the matrix CC^/A 2 . 
Thus, the entanglement entropy is 

S A = -tv{ PA \og{ PA )) = -(1 - A 2 £ aA log(l - A 2 £ a % ) - £ A 2 a 4 log(A 2 a,) 

= -A 2 log(A 2 ) £ a t + A 2 £ fli (l - log a,) + 0(A 3 ) . (18) 

Now, an explicit expression for the C matrix using standard perturbation theory is 

c nN = x (">W*bI°>°1 + om . (19) 

Eq + E — E n — E N 

Using this, the leading term in the entanglement entropy for small A is explicitly 

n^0,/V^0 C^O + ^0 - K - E N ) 

Interestingly, the entanglement entropy is not analytic in A at A = 0. Also, the leading 
order perturbative result (up to order A 2 terms which are not written explicitly) depends 
only on matrix elements of the interaction Hamiltonian between the vacuum and states 
with both subsystems excited. This follows since (|15p can be written as 

= (|0) + ^A n |n))®(|0) + ^^|0,iV))+ (C n ,N-A n B N )\n)®\N) . (21) 

In this expression, the entanglement would be zero without the second term, and in this 
term, C n N starts at order A while A n BN starts at order A 2 . The A and B coefficients 
do appear in the order A 3 contributions to the entanglement entropy (see Appendix [A}. 

Mutual information 

By a similar calculation, starting from a pure state in a theory with TL = 'H a (8)'Hb®Hci 
the leading contribution to I(A, B) in perturbation theory is 

\(N A ,N B ,N c \H int \0 : 0)\ 2 



(E -E Ni y 



/(A, B) = -A 2 log(A 2 ) 2 E + E 

I N A ^0,N B ^0,N c =0 N A ^0,N B ^0,N c ^0, 

(22) 

Similarly, when H = H A ®1-Lb® He ® Hd, at leading order in perturbation theory, 
the tripartite information I(A, B, C) is 

I(A, B, C) = +A 2 log(A 2 ) Y \(Na,N b NcNnlH^O)? + ^ ^ 

While I(A, B, C) can in general be positive, negative or zero, we see that the leading 
perturbative result for I(A, B, C) is always less than or equal to zero. This implies 
that to leading order in perturbation theory, 

I{A U B, C) < I(A, C) + I(B, C) . (24) 
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This result is not true for general systems^] Note also that if the matrix element of the 
interaction Hamiltonian between the free vacuum and states with all four subsystems 
excited is zerc0 then we will have 

I {A U B, C) = I (A, C) + I{B, C) . (25) 

to leading order in perturbation theory. In this case, the leading order contribution 
to mutual information between any two subsystems is completely determined from the 
mutual information between any pair of minimal subsystems^ 



3.2 Entanglement observables in quantum field theory 

For all of the observables discussed above, the essential quantities we have to compute 
are the reduced density matrices of the various subsystems. Given these quantities, 
we can compute the associated von Neumann entropies and mutual informations. In 
local quantum field theory, recent discussions of entanglement have focused on the den- 
sity matrices associated with bounded spatial regions. These are well-defined because 
(by locality) there are independent degrees of freedom in disjoint spatial domains, so 
the Hilbert space factorizes as required. The associated spatial entanglement entropy 
is typically divergent, even in free field theory, because in the continuum limit any 
spatial region contains an infinite number of degrees of freedom at arbitrarily short 
wavelengths. These divergences require regularization (e.g. by including an ultraviolet 
cutoff) and some care is needed to extract finite regularization-independent data [5]. 

Now, as discussed above, it is often more natural in quantum field theory to orga- 
nize degrees of freedom by momentum (or wavelength). Corresponding to any bounded 
subset of momenta in a field theory there are a finite number of degrees of freedom per 
unit spatial volume^] As a result, the entanglement entropy associated with such a sub- 
set should be finite for a finite volume system, increasing with the volume considered. 
For a translation-invariant state, we expect that the momentum space entanglement 
entropy will be an extensive quantity with a finite density S/V. We will verify this 
below. 

What observables can we compute? We can define the entanglement entropy S(P) 
associated to any subset P of the allowed moment at^l. the mutual information 7(P, Q) 
between any two subsets of momenta, or the tripartite entanglement I(P, Q, R) for any 
three subsets of momenta. We will focus on 

6 In particular, if A and B are completely uncorrelated, pab = PA® Pb, the opposite inequality, 
I (A U B, C) > I(A, C) + I{B, C) follows from strong subadditivity of entanglement entropy. 

7 In field theory, this will be true for theories with only cubic interaction terms. 

8 In field theory, such minimal subsystems will be the Hilbert spaces associated with modes at a 
single momentum. 

9 For a field theory at finite volume there will be a finite number of degrees of freedom in a bounded 
range of momenta. In the infinite volume limit, the set of allowed momenta becomes continuous. 
While there are now an infinite number of degrees of freedom with momenta in a finite region of 
momentum space, the number per unit spatial volume remains finite. 

10 More generally, P could represent a subset of the allowed single particle states. 
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• S(n), the entanglement entropy between all degrees of freedom with momenta 
above and below the scale \i. 

• S([[j,i, [12]), the entanglement entropy for a shell of momenta \i\ < \p\ < 

• S(p), the entanglement entropy for a single mode with momentum p. 

• I({\p\ < /Ui}, {\p\ > /U2}), the mutual information between degrees of freedom 
with momenta below a scale fii and degrees of freedom above the scale /i2- 

• I(p, q), the mutual information between modes with momenta p and q. 

These quantities probe the strength and extent of entanglement in momentum space. 

In free field theory, the Hamiltonian does not couple degrees of freedom with differ- 
ent momenta and thus all these measures of entanglement in momentum space should 
vanish. Adding a weak interaction term that couples degrees of freedom with different 
momenta modifies the ground state and should introduce a small amount of entan- 
glement between the various field modes. We can characterize this entanglement in 
perturbative quantum field theory by adapting the general results derived above. For 
the calculation of entanglement entropy, the two subsystems correspond to two com- 
plementary subsets A and A of the allowed momenta. The eigenstates \n, N) of the 
unperturbed Hamiltonian are elements of the Fock space basis |{nj}, {iV/}), where ri; 
and Nj are occupation numbers for particle states in the two subsets. The interaction 
Hamiltonian takes the form 



for some local Hamiltonian density "H/(x) that is polynomial in the fields and their 
derivatives. The matrix elements 



may be computed by expanding the interaction Hamiltonian density in terms of cre- 
ation and annihilation operators. The sum in ( 1201) is now over all states with at least 
one particle having momentum in the subset A and at least one particle having mo- 
mentum in the subset A. The nonzero matrix elements ( 126]) in the sum involve states 
with at most k momenta, where k is the number of fields in the interaction, and the 
momenta must add to zero since translation-invariance of the interaction Hamiltonian 
leads to a momentum-conserving delta function. 

The leading contribution ( 120]) to the entanglement entropy can be rewritten in terms 
of a projected two-point correlator of the interaction Hamiltonian (see Appendix). 
Below we will work directly with the expression ( 120]) . 




(K},{iv 7 }|tf 7 |o,o> 



(26) 
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4 Scalar field theory: entanglement between scales 



To develop some concrete examples of the general formalism, we will calculate momen- 
tum space entanglement entropy in d + 1 dimensional scalar theories with action: 

S = J d d+1 x{\{d,<j>f - \m 2 <j> 2 - ±r) ■ (27) 

For ease of formulation, we will first take the theory to be defined in a box of size L 
with periodic boundary conditions, and assume a UV cutoff at a scale A. 

We will compute the entanglement entropy S(fi) of degrees of freedom with mo- 
menta \p\ < fi with the high-momentum modes. Denoting by Pi and Pi the allowed mo- 
menta below and above /x, the Fock space basis elements are written as \ {n Pi })® {{np^). 
To use the general formula (T20I) for the leading contribution to the entanglement en- 
tropy, we need matrix elements of the interaction Hamiltonian between the free vacuum 
and the states with both low and high momenta excited. Recall that the fields can be 
expanded in terms of creation and annihilation operators as 

= (a p e-^ x + ale^) (28) 




where u p = \f p 2 + m 2 . The nonzero matrix elements for the interaction Hamiltonian 
between the Fock space vacuum and states with n particles excited^ are 



(Pi " " -Pn\HM = 1 5 f== (29) 



From ( 1201) . we then have 

SV> = -A> log(A*) £ + + 0(*) (30) 

{pi}n 

where the sum is over distinct sets of spatial momenta such that at least one momentum 
is below the scale /x and at least one momentum is above the scale \i. More generally, 
the entanglement entropy for the field modes with momenta in some set A is given by 
the same formula with the sum over distinct sets of momenta such that at least one 
momentum is in a set A and at least one momentum is in A. 

It is straightforward to take the limit of infinite volume. By the usual replacements 



E 



11 We are only interested in matrix elements between the vacuum and states with at least one low- 
momentum particle and at least one high-momentum particle. For tfi 3 and </> 4 field theory, the only 
such non-zero matrix elements have 3 and 4 particles excited respectively. For <j) n theory with n > 4, 
matrix elements with n — 2k > 3 particle states can also contribute, but for these theories we must 
also include <jf L ~ 2k counterterms in the action. For now, we focus on the case of (f> 3 and 4> 4 theory. 
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we find that the entanglement entropy density S(p)/V has a well-defined limit: 

S(p)/L d — — A 2 log(A 2 ) 7 J u [ T\d d Pl <5(pi / +, " +Pn) „ +0(A 2 ). 

(31) 

Here, the integral is again over distinct sets of momenta such that at least one momen- 
tum is below the scale p and at least one momentum is above the scale p. 

In practice, it is often simplest to calculate the derivative dS/dp, since the p- 
dependence comes only in the domain of integration, and this domain for S(p + dp) is 
almost the same as for S(p). In the difference 

S(p + dp) — S(p) 

the only contributions that don't cancel between the two terms are a positive contri- 
bution in which one momentum is in the range [p, p + dp] and all the other momenta 
have magnitude larger than p, and a negative contribution in which one momentum is 
in the range [p, p + dp] and all the other momenta have magnitude smaller than p. 

4.1 The 3 theory in 1+1 dimensions 

The simplest example is the (J) 3 theory in 1+1 dimensions J* 2 ! From (l3~Tj) . 



327H Jr x - LJ - LOi U1 2 UJ3{^1 + W 2 + UJ 3 ) 



Letting 

we find that 

1 dl 
2dp 



J(Pi,P2,p 3 ) = - f ■ ■ ^ (32) 

CU 1 CJ2UJ 3 {U 1 + UJ 2 + UJ 3 y 



poo rO 

j dpj(p,p,-p-p)- dpJ(p,p,-p-p). (33) 

J fl J-fi/2 



Evaluating the right hand side analytically for large and small p, we find tha 

v f.M/Ai _ (34) 



12/i 



{i + ln (^)} 



As discussed further below, the linear behaviour for small p is related to the linear 
growth in the number of degrees of freedom below scale p, while the falloff at large p 
is related to fact that a 3 is relevant in 1+1 dimensions so that the physics at large 
scales approaches that of the free theory, for which there is no entanglement between 
modes at different momenta. 



12 We work with a massive <fi 3 theory, so the theory is pcrturbatively stable. We can assume higher 
order interaction terms cj) 2n which stabilize the theory non-perturbatively but do not affect our leading- 
order perturbative calculations. 

13 To find we evaluate the expression for dl/d/j, and then integrate with respect to fi. The 
constant of integration is fixed by requiring that the entanglement entropy vanish as fx — »• 0. 
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Order A 2 terms 



Above we computed the 0(A 2 log(A 2 )) term in the entanglement entropy that dominates 
at infinitesimal A. At small, but finite A, the 0(A 2 ) term in (fT8|) could compete with 
this. To calculate this term we must determine the eigenvalues (and not just the trace) 
of the matrix CC^/A 2 , where 



C {Pi},{Pi} = "TTTTT^ , ^ P 't r^c i-r ( 35 ) 



and the sets {p^} and {Pi} must have either one and two elements or two and one 
elements. Thus the matrix M = C0 /A 2 has nonzero elements of the form M PA and 

M {puP2},{qim}- We haVe 

Mn „ = 5 pq —Y £±£±2 -. (36) 

' 8L ^ u p u Q u P {u p + ujq + u P y 



l P,Q 



Thus, for each p with \p\ < \i we have one eigenvalue 



*p = K T T^ 9 - (37) 

81 \P\>7lQ\>» ^p"q{"p + "P + ujqY 



The remaining block of the matrix M has entries 

5* 

j^j- , _ Pi+P2,qi+Q2 



Z {P1,J>2},{<?1,<?2} 



SL UJp 1 UJp 2 UJ qi UJ q2 L0p 1+p2 {UJ pi + U p , 2 + U pi+P2 ){u qi + U q2 + U pi+p2 ) 



where 5* indicates that we must have \p\ +P2I > H for a nonzero result. To find the 
remaining eigenvalues, we put M in block diagonal form, with one block for each P 
with \P\ > yU, where pi + P2 = Qi + Q2 = P- For the block labeled by P, we can label 
the matrix entries by p\ and qi, with 

w _ j_ 1 ngiMgi 



8L v /w Pl up_ Pl u gi up_ 91 (w pi + wp_ Pl + wp)(u;g 1 + wp_ gi + up) 8L 
where 

^(P) = ,. r • (38) 

^/0J p Up- p UJp{u p + Up- p + Up) 

A matrix of this form has only one nonzero eigenvalue, equal to 

*p = \ Y u7- (39) 

8L ^ U p UqUp{U p + U q + Up) 1 

We have have one such eigenvalue for each P with \P\ > p. 
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Figure 1: Leading contributions to S(/a) for 3 theory in 1+1 dimensions. Full result 
for S(fi) is proportional to A 2 (log(l/A 2 ) + 1) times bottom function plus A 2 times top 
function. 

Having found all the eigenvalues of CC^/A 2 , we can use ( fl8l) to write an expression 
for S(fi) including the order A 2 term. Recall that S(fi) = A 2 (l — log(A 2 ))^<2j — 
A 2 a i l°g( a i)- Taking L — > oo, 

E°*/ L = / tt I{pi) (40) 

and 

<k M^/L = J ^I( Pl ) \og(I( Pl )) (41) 

where 

fdp 2 dp 3 5{p l +p 2 +p 3 ) fAn . 
1 (PO = / ~Ta 7 , , \2 ■ 42 

Here, the asterisk indicates that p 2 < p 3 , and that p 2 and p% must have magnitude 
less than /i if p\ has magnitude greater than /i, while p 2 and p 3 must have magnitude 
greater than /i if has magnitude less than fi. 

We have plotted the two leading contributions (pLTj) and ( 140|) in figure [TJ We see 
that the two terms have a qualitatively similar behavior. In detail, the term ( )4T1) 
falls off slightly more slowly for large /i, behaving as 1/yU 3 log 2 (fj, 2 / m 2 ) compared with 
l//i 3 log(/i 2 /m 2 ) for ( T4*0l) . Thus, for fixed A and sufficiently large [i (of order m/X), 
the 0(A 2 ) term will be larger than the (9(A 2 log(A 2 )) term, although the qualitative 
behavior is similar. In this work, our focus is on the physics in the limit of small A, 
so in the remainder of the discussion we will concentrate (9(A 2 log(A 2 )) terms which 
dominate as long as we stay below the parametrically large energies of order 1/A relative 
to the mass. 
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(B) 



Figure 2: (A) Integration regions for <p 3 theory in 2+1 dimensions. (B) The function 
F(x) appearing in the entanglement entropy for 4 theory in 1 + 1 dimensions. 



4.2 The 3 theory in higher dimensions 

In general dimensions, the entanglement entropy for the modes below scale \x in the </> 3 
field theory is given by 

S(»)/L« = -A 2 log(A 2 )— * [ d d Pl d d P2d d Ps 6{P1+P 2 + P l ] . 2 +0(A 2 ) 

8{2ir) M J {p . }ft U^O^l^l + u 2 + cu 3 ) 2 

- - A2 MA 2 )^p^). (43) 



It is more convenient to compute 

1 dl d 
ujd-i^- 1 dfj, J B 



)dV((ftO),P,-(ftO)-p1 (44) 

A 



where J is defined in ( l32l) . and A and .B are the regions shown in Fig. [2j\ (symmetric 
between the vertical axis shown and the directions not depicted in the case d > 2). 
Here cod = 27i^ d+1 ^ 2 /T((d + l)/2) is the volume of the unit d-sphere. 
Explicitly, we have 

-J dp x j dp T UJ d -2PT~ 2 J(Px,PT) 

/(I poo 
dpx / / dp T u d - 2 p^ 1 J(Px,Pt) 

+ / dp x / dp T u d _ 2 PT~ 2 J{Px,PT) (45) 

J a JO 



u d ~i^ d ~ l d\i 



14 



where 

fi 2 + m 2 ^/p 2 x +PT + m 2 y/([i + p x ) 2 + p 2 - + m 2 

1 



m 2 



yj '[i 2 + m 2 + + p| + m 2 + \/{fi + Px) 2 +Pt + 
We find that in 2 + 1 dimensions, the entanglement entropy decreases with /i as 

27T 

JaOi) -> ^ (46) 
3/i 

when fi ^> m, while in 3+1 dimensions, we have 

/3(/i)^87r 2 (l-iln(2))/i 

for fi ^> m. We interpret the the 3 + 1 dimensional result as saying that in this 
case the /i 3 growth in the number of degrees of freedom below scale /x overwhelms the 
l//i falloff of the effective dimensionless coupling. These expressions are exact (and 
finite) as m — > 0. For 4+1 dimensions and higher, ( 1431) diverges - we will discuss this 
divergence below. 

4.3 4 theory 

Finally, consider the 4 field theory in 1+1 dimensions. From (]3~T]) . 



= -A 2 log(A 2 )— * / ggl±^±^j . 2 +0(A 2 - 
16(27r) J{p i \ IJL wi- -^(wi H hw 4 ) 2 



Thus, we study /(//) = f {p . }fi dp 1 dp 2 dp 3 dp i 5(p 1 + p 2 + p 3 + Pa) J(Pi,P2,P3,Pi), where 

<^(Pi)P2,P3,P4) _1 = Wi^^^^i + c^2 + ^3 + k^) 2 - It is again more convenient to 
evaluate 

ldi f f°° f pi f°° 1 

2 J- = W J dp 2 + J dpi J ^ dp 2 > J(pi,p2,n, -Pi - Pa - A*) 

- / dpi / dp 2 J(p 1 ,p 2 ,fi,-p 1 -p 2 -fi) = —F(n/m). 
A numerical integration determines F, giving the final result 

SM/V = -A 2 log(A 2 )^,^ / dxF{x) . (47) 
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The function F(x) is plotted in Fig. 03. By analyzing (analytically) the behavior of F 
for large and small x, we find that the entropy S(fi) behaves as /i/m 4 for small fi and 

as 



for large /i. As for the 3 theory, the decay at large \i is related to the fact that the 4 
theory in 1+1 dimensions is free in the UV. 

The leading perturbative contribution to the entanglement entropy S(n) of <f) A the- 
ory can be similarly evaluated in 2+1 dimensions. The integrals there are more dif- 
ficult to evaluate numerically, but are convergent. For 3+1 and higher dimensions, 
the integral expression for the leading contribution to S(fi) in the <f) A theory has a UV 
divergence, which we discuss further below. 

4.4 General remarks 

Massless limits: We found above that in two and higher space dimensions, the 
entanglement entropy S(fi) has a finite limit as we take the mass to zero. However, in 
1+1 dimensions, the results for both </> 3 theory and </> 4 theory diverge in the massless 
limit. These divergences suggest that S(fi) is not an "infrared-safe" quantity for a 
massless scalar theory in 1+1 dimensions. However, the ratio S(fi)/S(fio) has a finite 
limit if we hold /i and /xq fixed as we take m to zero. The result is 



Thus, while it may not be sensible to talk about S(fi) directly for m = and infinite 
volume, the ratio for different scales appears to be well-defined even in 1 + 1 dimensions. 

General understanding of large /i behavior: The results above agree with the 
following heuristic derivation of the power law behavior of S(fi) for large /i. The behav- 
ior is influenced by two significant effects. First, the number of degrees of freedom per 
unit volume below a momentum scale /i grows like fi d . All else being equal, we expect 
5* to scale like the number of degrees of freedom (for example, it is extensive). However 
the interactions in a general field theory depend on the scale, and this scale dependence 
also contributes to the behavior of S(fi). The dimensionless effective coupling for a <p n 
interaction at scale /x behaves as i / // i d + 1 - n ( rf - 1 )/ 2 . Since S goes like A 2 (plus logarithmic 
corrections), we can estimate that S(fi) should behave as 



up to possible logarithmic corrections. This is consistent with our results. At a techni- 
cal level this scaling arises in the integrals for entanglement entropy from two sources: 
(a) the measure factors (i.e. the density of modes), and (b) the energy denominators in 
the interaction terms. These are the same ingredients that affect the scaling of physical 
observables during renormalization. 



S ~ — In 2 (/i/m) 



5(//)/5(/x ) = (/V/i) 3 



(48) 
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Divergences: In various specific case considered above, we found that the leading 
perturbative contribution to S(n)/V is finite in the limit where the IR and UV cutoffs 
are removed. However, for the 3 theory in 4+1 and higher dimensions or the </> 4 the- 
ory in 3+1 and higher dimensions, the integral expressions for the leading perturbative 
contribution to S(/j) diverge. The divergence is associated with the sum over states 
in (KOj) . or in the sum or integral over the momenta in (130 j) or (l3"Tj) respectively. The 
leading divergence comes from the sum over states where one momentum has magni- 
tude less than /j while the rest have magnitudes greater than \i. The divergence is 
proportional to a power (or logarithm) of the UV cutoff A. 

Of course, ultraviolet divergences are commonplace in quantum field theory. Typi- 
cally, they are associated with integrals over momenta that appear beyond the leading 
order in perturbative calculations, and are dealt with by expressing the results in terms 
of renormalized (physical) parameters rather than bare parameters. However, here the 
divergences appear in leading order perturbative results. Since the bare and renormal- 
ized parameters are the same at leading order in perturbation theory, the divergences 
will not be eliminated by expressing the results in terms of renormalized parameters!^! 

Such divergences in leading order expressions indicate a breakdown of perturbation 
theory for the specific quantity that is diverging. To see this, note first that similar 
divergences appear even in fermionic theories, for example fermions in 1+1 dimensions 
with a (ipijj) 2 interaction (see Appendix C for details). Furthermore, the divergence is 
present even at finite volume, since it is associated with the infinite number of high- 
momentum modes which are still present with an IR regulator. But for a theory of 
fermionic fields at finite volume, the Hilbert space associated with degrees of freedom 
with momentum below a scale /x is finite-dimensional. In this case, there is an upper 
bound S(fi) < \og(N), where N is the dimension of the Hilbert space. Now, consider 
the theory with a UV cutoff A. The leading perturbative expression for S(/a) will be 
finite for any finite A. But since this expression diverges as A is taken to infinity, there 
will be some finite A above which this leading contribution to S(fi) is larger than the 
bound \og(N). Here, A is still finite, so S(fi) is clearly well-defined, and the correct 
result for S(fi) must certainly be less than log(iV), so the only possibility is that the 
leading perturbative expression is not a good approximation to the correct result. 

Our conclusion should not be particulary surprising: regardless of how small the 
coupling parameter of a theory is, there will always be quantities that cannot be com- 
puted in perturbation theory. Here, the breakdown of perturbation theory seems to 
be associated with computing the entanglement entropy between a finite set of modes 
with the infinite set of degrees of freedom above scale \x. We will see below that in 
cases where perturbation theory breaks down for this quantity, it is still possible to 
perturbatively calculate less inclusive quantities, such as the mutual information be- 
tween degrees of freedom associated with two finite regions of momentum space. In 
cases where no divergence appears at leading order, the finite leading-order perturba- 

14 We do expect the standard divergences to appear in higher order perturbative corrections, even 
for quantities whose leading order result is finite. These divergences should be cured in the usual way 
by expressing results in terms of physical parameters, or by using renormalized perturbation theory 
with the appropriate countcrterms. 
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tive result should be reliable so long as subsequent terms in the perturbative expansion 
(after renormalization) are small compared to the leading terms. 



5 The extent of entanglement between scales 

So far, we have considered the entanglement between modes in a field theory above 
and below some scale /x. In this section, we ask about the entanglement entropy 
associated with a single mode of the field theory, or the mutual information between two 
individual modes. A version of the former observable has been considered previously in 
the condensed matter literature (see e.g. [7]). We can also consider the entanglement 
entropies of bounded regions of momentum space. These sorts of observables are useful 
for two reasons: (a) they can be finite even when the entanglement entropy for the low- 
energy density matrix diverges, (b) they are a much more sensitive and clear probe of 
the extent of entanglement since they don't sum over the entire tower of UV modes. 

5.1 An aggregate measure of the range of entanglement 

The quantity S(p,) measures entanglement between the complete set of degrees of 
freedom below the scale \i and the complete set of degrees of freedom above the scale 
\i. Is this entanglement largely between modes just above and below the scale /i, or is 
the entanglement "long-range" in momentum space? 

One way to address this question is to consider the entanglement entropy for the 
annular region \i\ < \p\ < \ii in momentum space. If the entanglement is short-range, 
then for /i 2 S> /Ji, the entanglement entropy 5([/ii, /i 2 ]) = S{p,\ < \p\ < fi 2 ) should be 
dominated by entanglement between modes just above and below the scales /ii and /i 2 - 
In addition, these separate contributions to the entanglement entropy should be well 
measured by S^/^) and S(fii). Thus, for short-range entanglement we would expect 



Alternatively, consider the mutual information between the degrees of freedom with 
\p\ > ^2 and \p\ < jii: 



For short-range momentum space entanglement we expect (149]) . Hence, when /i 2 3> /ii 
we expect that /(/ii,/x 2 ) ~ 0. The rate of falloff of I(/i 2 , Hi) as ii 2 //ii increases from 1 
is a characterization of the extent of entanglement. 

In <f) 4 theory the infinite volume expression for 7(/ii,/i 2 ) is (using ( 1221) ) 



A*2 > A*i • 



(49) 



J(/ii,/i 2 ) = S(ni) + S(n 2 ) - Sfl/ii,/^]) • 



(50) 



S([^^})/V = -\Hog(\ 2 )^ [H 



d d pi (2vr) d 5(pi +P2+P3 + Pa) 



+ 0{\ 2 ) 



2(27T) d {Ux +U2+UJ3 + UJa) 2 ^1^2^3^4 



(51) 
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where the asterisk indicates that we integrate over momenta such that at least one 
\p\ is in the range [pi, p 2 ] and at least one \p\ is not in this range. For simplicity, we 
specialize to d = 1 and take the mass m = 1. It is simplest to first evaluate the quantity 

* S (52) 



d\i\dp 2 



We can see that the only contribution to this will be from regions of the integral above 
where one momentum is at p\ and another momentum is at p 2 . 

This is equal to the integral above with p\ = ±p 2 , p 2 = an d the remaining |p|s 
either both inside or both outside the interval [pi, p 2 \. The distinct choices of momenta 
satisfying these constraints are 

Pi = P2=Pi P3 6 (-oo, -2/i 2 - Pi] U [-2/ii - p 2 , + ^2)] 

Pi =1^2 P2 = -^l P3 6 (-00,-2^2 +A*l] U [-Atl,|(/il -/i 2 )] 



or momenta obtained from these via pj — > — p^, where in all cases, P4 is determined by 
the delta function constraint. Thus, we have 

-A 2 log(A 2 )^ ifi ^ { f'^ 2 '^ 1 dp J(p 2 , p h p, -p - p! -p 2 ) 



V dp x dp 2 y 1216(27r) 



+ /^-/tf 2) dp J p> -p- hi- ^2) 

+ /rI M2+m <^ ^2, -p+yi- M2) 
+ J-lT"^ dp J ^ 2 ' -p+ hi- ^2)} 



where 



■J(Pl,P2,P3,Pi) 



1 



+ Cd 2 + ^3 + W4) 2 WlW2^3^4 

To determine S(p±, p 2 ) from this expression, we can use S(p, p) = 0, S(0, p) = S(p), 
[ J^(0,/z) = jt-(p). From these, we have 

dS , , d 2 S , dS" 



<9/J : 



/■mi d 2 S dS 
-(^1,^2)=/ dp l ———(p u p 2 ) + —(p 2 ) (53) 
J dp x dp 2 dp 



and 



n r 1 * ds , _ /■« d 2 s .„ „ . p» ds,„. 

b{pi,p 2 )= / 07/23 — (^1,^2) = / "^2 / a/ii- — - — {p 1 ,p 2 ) + / dp 2 — {p 2 ). 
Jm dp 2 Jo dp x dp 2 y w d/x 

(54) 

Here, S^/x) is the quantity that we evaluated in previous sections. 

To investigate whether the entropy iS(/xi, /x 2 ) is dominated by entanglement between 
degrees of freedom close to \i\ and p 2 , we can vary p 2 and ask whether the variation 
of S(pi,p 2 ) is well approximated by the variation of S(p 2 ) (these variations would be 
equal if S(pi,p 2 ) = S(pi) + S(p 2 )). From (J53"]) . the difference between the variations 
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(A) 



Figure 3: Ratio of first and second terms in ([331 vs /i = (fi 2 — for (A) fi\ = 1 

and (B) fii = 4. This is a measure of the range of entanglement in 4 theory in 1 + 1 
dimensions. We have taken the mass to be m = 1. 

is equal to the first term on the right hand side, so we ask whether this term is small 
compared with the other term. In Fig. [31 we plot the ratio of these terms as a function 
of /i = (/i 2 — yUi)/m. The ration declines as 1/ln/x increases and approaches a finite 
value as (/i 2 — Hi)/m — > 0. 

The slow rate of decline is surprising given that the 4 theory in 1+1 dimensions 
enjoys the property of decoupling. Note however, that the quantity we are computing 
integrates over all of the UV modes. Thus it is an aggregate measure of entanglement. 
A more refined way to ask about the range of entanglement in momentum space is to 
consider the mutual information between individual modes at two different momenta 
p and q as we do below. We will see that this mutual information falls off as a power 
law with \q\ when \q\ 3> |p|. 

5.2 Single mode entanglement 

In this section, we calculate the entanglement entropy for a single mode with momen- 
tum p. This measures the entanglement between a single mode and the rest of the field 
theory. The leading result for a <fi n scalar field theory follows immediately from ( 130]) : 



s(p) = -A 2 io g (A 2 ) y 9 , ^2±z+ 

{P2,.",P™} 




where the sum is over all distinct sets of (n — 1) momenta^ In the infinite volume 
limit, this gives 



S(p) = -A 2 log(A 2 ) , * [ f]d d Pi 5d(j) 




15 For p = 0, we have the further restriction that not all momenta are zero. 
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By rotational invariance, the result is a function only of \p\. All explicit factors of the 
volume have canceled out without dividing by volume on the left side. 

A natural interpretation of this finite quantity is that it gives the entanglement 
entropy density for degrees of freedom in an infinitesimal range d d p around the mo- 
mentum p. The number of modes in the box d p is proportional to spatial volume, 
so if the entanglement entropy for one mode has no explicit volume dependence, the 
entanglement entropy for the set of modes in the box should be proportional to vol- 
ume. This entropy is also proportional to the momentum space volume d d p (if this 
is infinitesimal), so the entanglement entropy associated with an infinitesimal volume 
d d x in position space and volume d d p in momentum space takes the form 

dS{p) = s i(bD • ( 56 ) 



It is interesting that the phase space volume appears naturally hereJ 16 ! 

As an explicit example, we have plotted Si(p) for 3 theory in two, three, and four 
spacetime dimensions in Fig. HJ In the figure, the entropies are normalized by their 
value at p = 0. For 1+1, 2+1, and 3+1 dimensions, the entanglement entropy decreases 
like l/p 4 ,l/p 3 , and 1/p 2 respectively. Thus we see that in this case the entanglement 
of a single mode with the rest of the theory declines as power-law of the momentum, 
even though we found above that the integrated entanglement between modes above 
and below scales fi2 and /ii only declines logarithmically. The slow decay in the latter 
case is arising from the sum over modes. 

5.3 Mutual information between individual modes 

It is also interesting to investigate the mutual information between two specific field 
theory momenta. In the large volume limit, the natural quantity to consider is the 
mutual information between degrees of freedom in an infinitesimal range d d p around 
some momentum p and degrees of freedom in an infinitesimal range d d q around some 
momentum q. Starting from the basic formula ( )22|) . only contributions from the second 
term in curly brackets survive the large volume limit. These involve matrix elements 
between the vacuum state and states where one particle is excited in each of the regions 
d d p and d d q, and the remaining particles lie outside these regions. The resulting mutual 
information is proportional to d d p and d d q and to spatial volume, so we have: 

I(p,q)/V = d d pd d qX(p,q);. (57) 

For <p n scalar field theory in d + 1 dimensions, the leading contribution to X is 

T(p,q) = -A log A ) / I d pi ■ ■ V2+0(A 

2 n (2Tx) d ^ i) J {p3> ... >Pn} f = \ uj p u q u 3 ■ ■ ■ u n {u p + ■■■ + u n ) 2 

(58) 
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Note that while this entropy is spatially extensive, it is not extensive in momentum space. That 

r' 



it, it is not true that S(R)/V = J R d d pf(p) 
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2 4 6 8 10 



Figure 4: Single- mode entanglement entropy vs magnitude of mode momentum for 3 
field theory in 1+1 (bottom), 2+1 (middle), and 3+1 (top) dimensions. The entropies 
are normalized by their values at p — 0. 



where the integral is over distinct sets of n — 2 momenta. For 3 theory, this is 

Z(p,$ = -A 2 log(A 2 ) , * ^ + 0(X 2 ) (59) 

Thus, the mutual information is enhanced when p, q or (p + q) are near zero, and for 
fixed p, the mutual information falls off as 1/|<?| 4 for large q. While this expression 
gives the formal leading order result in any number of dimensions, we will see below 
that it should only be trusted as an accurate approximation to the exact result for 
d < 4 space dimensions. 



5.4 Convergence and validity of leading order expressions. 

As for the entanglement entropy S(fi) considered in the previous section, the integrals in 
the leading order contributions to the mutual information and entanglement entropy of 
single modes can contain UV divergences. As we argued in Sec. I4.4[ such divergences 
indicate a breakdown of perturbation theory for the quantity in question. In this 
subsection, we classify the scalar field theories for which the perturbative calculation 
of single-mode quantities si(p) and T(p,q) gives sensible results. 

We begin with the expression (158]) for the single mode mutual information of <ft n 
scalar field theory in d+ 1 spacetime dimensions. Naively, this will converge (i.e. there 
are enough powers of momenta being integrated over in the denominator) if 

d < 1 + . (60) 

n — 6 
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Thus, we have convergence in any dimension for n = 3 for d < 3 for n = 4, for d < 2 
for n = 5, and only for d = 1 for any higher n. 

Since higher order interactions (i.e. interactions with more powers of the field) are 
more likely to lead to divergences, we should be concerned that such higher order in- 
teractions generated in the quantum effective action will produce divergences at higher 
orders in perturbation theory. In 3 theory we get an effective n vertex at order 
A n from a one loop diagram. As a function of the external momenta, this scales like 
-^jp2n-d+i ag these momenta are taken large. The contribution to I(p,q) from such a 
vertex will naively be convergent if 

d<5 + —?— . (61) 
n — 1 

This is satisfied for any n so long as d < 5, but leads to a divergence in 6 and higher 
space dimensions. Thus, it appears that I(p, q) can be computed in perturbation theory 
for 3 theory in d < 5 (the same dimensions for which the theory is renormalizible) . 

For 4 theory, the effective action contains effective 4> 2n interactions coming from 
one-loop diagrams at order A n . These scale with external momenta like l/p 2n +1 . The 
contribution to I(p, q) from such a vertex will naively be convergent if 

d < 3 



In- 1 

which is satisfied for any n as long as d < 3. Thus, it appears that I(p, q) is a well 
defined quantity for 4 theory in d < 3 (again, the same dimensions for which the 
theory is renormalizible). 

An almost identical analysis shows that the mutual information between degrees 
of freedom in any two finite region of momentum space converges whenever I(p, q) 
converges. Note that it would be incorrect to suppose from the considerations above 
that the leading order I(p, q) is necessarily well defined for every renormalizable theory. 
For example, according to (I6ip . the leading order I(p, q) diverges for the renormalizible 
4> 6 theory in 3 dimensions. We can also ask when the entanglement of a single mode 
(or a finite region of momentum space) with the rest of the field theory is well defined. 
For <p n theory, we find convergence for 

, 2n - 1 

d < — (62) 
This result extends to entanglement entropy of any finite region of momentum space. 



6 Comments 

We have obtained a number of results for the scaling of entanglement entropies and 
mutual informations with the upper bound on a momentum interval. At a technical 
level, these results all follow from the density of modes (measure) in the integrals over 
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Theory 


Dimensions where I(A,B) converges 


Dimensions where 5(A) converges 




d < 6 


d < 4 


4 


d < 4 


d < 3 


5 


d < 3 


d = 2 


0n>6 


d = 2 


d = 2 



Table 1: Dimensions where momentum space mutual information and entanglement 
entropy converge. The results apply for any bounded regions A and B in momentum 
space. 



momenta, and the energy denominators in the interactions. These are the same ingre- 
dients that lead to the decoupling property of local quantum field theories. Indeed, 
decoupling is usually understood simply as the power law suppression of higher di- 
mension operators in a low energy effective theory. This suppression means that high 
momentum degrees of freedom have weak effects on the dynamics of low-momentum 
degrees of freedom other than renormalizing the interaction strengths and wavefunc- 
tions. Our study of entanglement between degrees of freedom with different momenta, 
and the resulting entanglement entropies and mutual informations, refines this under- 
standing of the influence between momentum scales. 

In more detail, "decoupling" between UV and IR physics implies that starting from 
a generic action S A (gj) at scale A, that depends on an infinite number of parameters 
gj, the Wilsonian effective action at a much lower scale /i will be very close to some 
action S^(gi) in a family parameterized by a small number of physical parameters 
gi. In other words, the operation of integrating out degrees of freedom to successively 
lower scales results in a flow in the space of Wilsonian effective actions that converges 
to a low dimensional subspace at scales /i<A. Now, according to ([5]), the Wilsonian 
effective action S!^ at scale /i completely determines the reduced density matrix 
for the degrees of freedom with \p\ < fi. Thus we conclude that for the ground state 
of a generic field theory defined at scale A, the reduced density matrix for the degrees 
of freedom below some much lower scale \i will be very close to some family of density 
matrices p(/i, g^) that depend on a small number of physical parameters gi. Conse- 
quently, knowing the state of the low-momemtum degrees of freedom tells us relatively 
little about the details of the state at much higher scales. 

The paucity of information about UV physics contained in the low-momentum 
density matrix should be reflected in some of the measures of quantum information 
we have discussed. Specifically, it seems likely that there is a connection between the 
decoupling behavior of field theories and the power-law fall off in mutual information 
observed in Sec. It would be interesting to make this connection precise. 
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Relation to AdS/CFT: In the context of gauge-theory / gravity duality (the 
AdS/CFT correspondence) [10], there is now evidence that certain measures of en- 
tanglement in quantum field theory carry geometrical information about the dual 
spacetime. For field theories with a weakly curved dual gravity description, Ryu and 
Takayanagi have proposed rt that the entanglement entropy for a spatial region A is 
proportional to the area of the minimal surface A in the bulk space whose boundary 
coincides with the boundary of A, 

_ Area(A) 

While the proposal has not yet been proven, it has passed a variety of checks (see e.g. 

m na eh us])- 

Given the holographic interpretation of position-space entanglement entropy, it is 
natural to ask whether the momentum-space quantities considered in this paper have 
some simple dual geometrical interpretation for field theories with gravity duals. As 
an example, the quantity S(fi) measures the entanglement between degrees of freedom 
above and below the scale /i. Since energy/momentum scale in holographic field theories 
corresponds to radial position in the dual geometry, we might guess that S(fi)/V is 
related to the area (per unit field theory volume) of a surface separating the IR region 
r < r(fi) of the dual geometry from the UV region r > r(/i). For the dual geometry to a 
translation-invariant field theory state, this area function is a well-defined observable0 
However, we currently have no way to check whether this or some similar observable 
corresponds to momentum-space entanglement entropy, since we cannot calculate this 
entropy for any strongly coupled field theory with a gravity duall^l 



Relation to DMRG and MERA: Here we have explored various aspects of entan- 
glement in quantum field theory and the connection to renormalization theory. In the 
condensed matter literature, the ideas of entanglement and renormalization have come 
together previously in various schemes for approximating the ground state of many- 
body systems [?.I19]. While the focus and details of that work are rather different from 
the present discussion, it may be useful to briefly review those ideas here. 

Consider a quantum many-body system described by some lattice of degrees of 
freedom, for which the Hilbert space decomposes as a tensor product of Hilbert spaces 
for the individual sites. The dimension of the full Hilbert space is d N where d is the 
dimension of the individual Hilbert spaces and N is the number of sites. A general 
state (and in particular, the exact ground state of the system for a given Hamiltonian) 

17 If the spatial part of the dual metric is dr 2 + f(r)dx 2 , the area of the surface at radius r(/u) per 
unit field theory volume is proportional to a power of /(r(/i)). 

ls S([i) will probably not always correspond in a simple way the specific area observable mentioned, 
since that area would be well-defined even for gravity duals of 0+1 dimensional field theories, for 
which there is no way to divide up the degrees of freedom by spatial momentum, and therefore no 
way to define S(fi). Of course such low-dimensional gauge/gravity dualities (e.g. AdS2/CFTi also 
have many other special features [T()l IT7] ). 
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can be represented exactly by a tensor T ll "' lN that gives the coefficient of the basis 
state ® ■ ■ ■ <g) |ijv)- 

A general numerical determination of the ground state is impractical due the the 
large number d N of independent coefficients. For certain systems, usually in 1+1 
dimensions, an efficient variational approach to approximating the ground state is to 
consider tensors T that can be decomposed into contractions of lower-rank tensors. For 
example, the "Matrix Product State" (MPS) decomposition corresponds to 

= (M 1 )i 1 ia2 (M 2 )jf 203 • ■ ■ (M N )Z ai . 

In practice, one uses this decomposition as a variational ansatz, varying the individual 
matrices M % to arrive at the best approximation to the ground state Cite. If the 
dimension of the matrices M % is large enough, any tensor T can be represented in 
this way, so the variational method gives an exact result. However, for a wide class 
of systems, it has been found that the ground state can be well approximated by 
matrices of much lower dimension. In this case, the matrix product ansatz represents 
a truncation of the Hilbert space to a subspace of lower dimension, and in cases where 
it is effective, the true ground state is close to the ground state in this subspace. 

It turns out that the success of this method is related to the entanglement properties 
of the ground state. The optimal method of truncating to a lower-dimensional Hilbert 
space is to retain as much of the entanglement entropy for the various subsystems 
(blocks of sites) as possible^ The procedure works most efficiently (i.e. for smallest 
matrices M) when there is limited entanglement between the subsystems corresponding 
to blocks of sites. For systems with a highly entangled ground state, the method is 
much less efficient. 

Another approach that is more successful in cases with long-range entanglement is 
the "Multiscale Entanglement Renormalization Ansatz" (MERA [19]). In this approach 
the tensor T is represented by an iterative procedure. The tensor is first written in 
terms of a "disentangled" tensor T using unitary matrices U : 

rpil—ilN _ (Tj( n )\ili2 /rr(n)\*2JV-l«2iV rfijl—j 2 N 

1 (n) - \ U 1 )hh ' ' ' \ U N )j 2 N-d2N 1 (n) 

and then T is represented in terms of a lower rank tensor using "projectors" P: 

rpjl—j2N _ ( p( n )\jlj2 ( r>(n)\j2N-lj2NrpI 1 --I N 

1 (n) — V-M )h "'\ r N )l N 1 (n+l) • 

The latter step can be understood as a "coarse-graining" of the system, though the 
dimension of the index space I is not necessarily the same as that of the original index 
space i. The original tensor is thus represented by the individual matrices {U^) 1 ^,-, 

and (P} Yj which are the variational parameters used to approximate the ground state. 

19 This idea arose first in the "Density Matrix Renormalization Group" (DMRG) [18] an iterative 
renormalization procedure on the state of the system that truncates the Hilbert space in each step 
while retaining as much entanglement entropy as possible. The DMRG is now understood to give 
results equivalent to this MPS variational method. 



26 



The introduction of the U matrices is motivated by the observation that coarse 
graining works most efficiently when there is little entanglement between the adjacent 
blocks. The unitary matrices U can remove short-range entanglement between adjacent 
blocks before coarse graining. In this way, the matrices Ui n ) encode the entanglement 
between sites at the nth level, which corresponds in the original picture to blocks of 
2 n sites. Thus, in the MERA representation of a ground state the unitary matrices U 
encode entanglement at different scales. This information is certainly related to the 
scale-dependent entanglement entropies considered in this paper, though the MERA 
entanglements would seem to be more closely related to position space entanglement. 
Also, the original MERA applies only to discrete systems, though an extension to 
continuum quantum field theories has been recently proposed in [20]. An interesting 
connection between MERA and the AdS/CFT proposal above has been given in [21]. 
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A Entanglement entropy at 0(A 3 ) 

The calculations in Sec. 13. II can be extended to give an expression for the (9(A 3 ) terms 
in the perturbative calculation of entanglement entropy for a general system. We find 
that after a similarity transformation, the density matrix f fl6|) can be written as 

PA ~{ act - ABC \ - CB^ti J +U ^)- ^3) 

The eigenvalues of this matrix (up to corrections of order A 4 ) include the top-left 
element of the matrix, and the eigenvalues of the lower-right matrix. At leading order, 
the lower-right matrix is C\C\ where C\ is the order A term in C. We defined the 
eigenvalues of this matrix to be aj. Finding the eigenvalues of the lower- right matrix 
after the higher order terms are added is a problem formally equivalent to ordinary 
time-independent perturbation theory in quantum mechanics, so we can express the 
result in terms of the eigenvalues and eigenvectors of C\C\. 

Using this approach, the result for the entanglement entropy up to order A 3 is 

SQi) = A 2 (-log(A 2 ) + l)tr(a 1 C 1 t )-A 2 ^a i log(a J ) 

i 

+A 3 (- log(A 2 ))tr(C 1 C 2 t + C 2 C\ - A 1 B 1 C\ - C X B\A[) 
-A^logCa^ldCt + C 2 C{ - A.B.Ct - CMAh) 

i 

where a« and u$ are the eigenvalues and eigenvectors of the matrix CiC} and the 
subscripts indicate the order in perturbation theory. 
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B Momentum-space entanglement and correlators 

Starting with the general expression ( l20l) for the leading order perturbative contribution 
to entanglement entropy, we can now specialize to the case of quantum field theory. 
We find that 



/oo 
rfrr(0, 0|#/|n, ^e^ - "^^^^, iV|#j|0, 0) + C(A 2 



/oo 
rfrr(0, 0\e HoT H ie - HoT \n, N)(n, N\H T \0, 0) + C(A 2 ) 

/■oo 

= -A 2 log(A 2 ) / c/rr(0,0|e^ or i/ / e-^ or n Aj fi" / |0,0) + C(A 2 ) 

/•oo 

= -A 2 log(A 2 ) / drT(^(-iT)n A F 7 (0)) + G(A 2 ) 

rfrr / d z xd i y{H I {-iT,x)Ii A 'H I ^,y)) + 0{\ 2 ) 

/OO /» 
rfrr J d 3 x(%(-zr,x)n^j(0,0)) + O(A 2 ) 

Here, we use the standard definition of time-dependent operators in the "interaction 
picture" : 

#j(t) = e iHot H ie - iHot . 

The operator II is projects to intermediate states with at least one particle having 
momentum in the subset P and at least one particle having momentum in the comple- 
mentary subset of momenta. 

The factor of volume in the last line comes from the y integral in the previous line, 
which is trivial since the correlator in that line can depend only on the combination 
x — y. The entropy per unit volume S(P)/V will have a finite limit, so that S(P) is 
an extensive quantity. 



C Entanglement entropy in a fermionic system 



Here we calculate the entanglement entropy in a fermionic theory with a (ipip) 2 inter- 
action. Consider for definiteness the renormalizable theory in 1+1 dimensions. The 
fermion fields are expanded as 



^0) = Y] tt^= 

- L-2 y/ZUp 



-ipx 



blv(p)e ipx ) 



(64) 



As a straightforward application of (120]) the entanglement entropy is 



^ = -A 2 log(A 2 )££ 



<{p,0i,-,{p,04 (^) 2 o) 



(oji + u 2 + u 3 + u A y 
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0(A 2 



(65) 



where t indicates the type of fermion (i.e. particle or antiparticle). The star indicates 
that the sum over momenta is restricted to the set where at least one momentum is 
above and at least one momentum is below the scale fi. Substituting the expansion 
(JMD into (ESD 

6 4 * S 

Sf, = -A 2 log(A 2 )-— V— — \u(p 1 )v(p 2 )u(p 3 )v(p 4 ) - u(p 1 )v(p 4 )u(p 3 )v(p 2 )\ 2 

(66) 

where 6 is the number of ways or choosing 2 particles and 2 antiparticles. Using 1+1 
dimensional spinor and gamma matrix identities, and passing to the infinite volume 
limit we are left with 



6 f* j_ j„ r/V^ „ \ G°i ' Vz - m 2 )(p 2 ■ Pa - m 2 ) 

T 



S^/L = -X 2 log(X 2 )^y 3 J* d Pl ...dp 4 5(£p 



In the region where three momenta are taken to be large this integral diverges linearly. 
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